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Abstract. We present a model in which the parameter w approaches — 1 near a particular value of z, 
and has significant negative values in a restricted range of z. For example, one can have w ~ — 1 near 
z = 1, and w > —0.2 from ^ = to z = 0.3, and for z>9. The ingredients of the model are neutral 
fermions (which may be neutrinos, neutralinos, etc) which are very weakly coupled to a light scalar 
field. 

This model emphasises the importance of the proposed studies of the properties of dark energy 
into the region z > 1 . 
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INTRODUCTION 

About twenty five years ago the possibility of neutrinos interacting by exchange of 

very light neutral scalars to produce interesting cosmological and astrophysical effects 
was considered by Kawasaki, Murayama and Yanagida [1], and Malaney, Starkman and 
Tremaine[2]. 

Three of us later investigated a similar system, in which the scalar particle has a mass 
of order 1/a.u., which leads to neutrino clustering [3]. 

Following the observation of the re- acceleration of the expansion of the Universe [4, 
5], which led to the concept of dark energy (for a recent review see ref. [6]). Fardon, 
Nelson and Weiner [7] noted that the inferred energy density of dark energy, of order 
(2.4 meV)4 [6], was, remarkably, of the order of the experimentally inferred value of 
neutrino masses. They also utilized the concept of a scalar field interacting with the 
neutrinos to obtain the desired effect. 

It is worth emphasizing that negative pressure is not a strange feature in physics, in 
that any self bound system has an equilibrium density, and will have a higher energy per 
particle as the density is decreased, and thus a negative pressure in this region. Because 
this feature exists in our neutrino clustering model, we realized that it could produce 
negative values of w and was thus relevant to the dark energy problem. We provided 
some early comments [8], and then showed how it led to a system with w ~ — 1 in 
a restricted range of the development of the universe [9], with w ^ 1/3 at very early 
times, and w — near the present. In this paper we will be elaborating our model of 
dark energy. 

To obtain the energy density of the neutral fermion - scalar field system as a function 
of z one needs to make assumptions about parameters. Typically w has a minimum as a 



Epoch Dependent Dark Energy 



August 10, 2009 



1 



function of z, and the position of the minimum may be moved by a choice of parameters, 
but a rapid variation of w between the minimum and the present is characteristic of our 
model. This lends support to attempts to probe the z dependence of dark energy in more 
detail. 

It is worth noting that Mota et al[10, 11] have recently introduced the possibility of 
neutrino clustering in the neutrino - scalar field model to obtain a dependence of the 
dark energy density on z. Our work shows that a variation of w and the dark energy 
density pE with epoch can be obtained in a homogeneous model. 

After outlining the model, we display our results for w as a function of the density 
of the neutral fermions, and then discuss how to choose the parameters of the model to 
obtain results for w as a function of z, and discuss the implications of the results. 



THE NEUTRAL FERMION - SCALAR FIELD MODEL 

Following ref. [3], the equation of motion for a neutral fermion field, ^^f, interacting with 
a scalar field, ^, is 



[d^ + m]](l> = g\\f\\f (1) 
\\f = -g(^\\f, (2) 



with h = c = 1 . The nonlinear scalar selfcouplings are omitted here, even though they 
are required to exist by field theoretic selfconsistency [8], as they may consistently be 

assumed to be sufficiently weak as to be totally irrelevant. The parameter m^^ is the 
renormalized vacuum mass of the isolated neutral fermion, and takes into account the 
contributions from the electroweak theory, as well as contributions from the vacuiun 
expectation value of the new scalar field 0. 

We look for solutions of these equations in infinite matter which are static and transla- 
tionally invariant. In the Freidman equation description of the evolution of the universe 
we use the adiabatic approximation, assuming that the scalar field expectation value 
takes the value determined by the fermion density at the appropriate time^ Equations 
(1,2) then lead to an effective mass for the neutral fermion of 

2 

m: = mf-^W- (3) 

These equations are simply the equations of Quantum Hadrodynamics [12]. 

As in [3], we act with these operator equations on state which is a filled Fermi 
sea of neutral fermions, with a fermion number density p„ and Fermi momentum kp, 
related, for Majorana particles (which we assume henceforth), by p„ = kp^ /{3k^). We 

introduce the parameter ^0 = / (;r^m^) , and the variables}; = m*Jm^n^ , 



We have studied the corrections to the adiabatic approximation in the fermion- scalar field model and 
found them to be small, but have not yet investigated the corrections to this approximation when that 
model is coupled to the Freidman equation. 
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X = k/m^\xF = kF/m^\ and ep = \Jx\+y^. Equation (3) becomes an equation 
fory: 



y=\- 



epXF-y In I J 

,(0) 



(4) 



The total energy of the system, E = em„ N, with N the number of fermions, and 



Koxi 



(5) 



As xf goes from to oo, j varies monotonically from 1 to 0, behaving as 
y Ki I — (^ox^) /3 for small xf, and y ~ 2/ [Koxjr) for large xf- 

For the fermion system to be bound, the minimum of e must be less than 1. We 
find that for sufficiently large the neutral fermion - scalar field system is bound [3], 
and the total energy displays the characteristic behaviour of a self-bound system noted 
above, suggesting that we should see negative pressures in this system. We now go on 
to compute the equation of state parameter. 



DETERMINING w 

The Equation of State, or the equation for the pressure, P, as a function of the other 
variables of the system, is represented through the equation of state parameter, w, which 
is defined by the proportionaUty between P and Pe, the energy density 

P = wpE. (6) 

As is well known, for radiation, w = j, for cold non-interacting matter, w = 0, and, for a 
Cosmological constant, w — — I. Any w < describes a system with negative pressure. 

Using the equations of the Friedman-Lemaitre-Robertson- Walker (FLRW) metric, 
with a scale parameter a, the equation of state parameter is shown to satisfy 

1 dlnpE 

l+w=--^- . (7) 

3 oma 

In these equations pg is the total energy density. Assuming for this paper that the only 
source of this energy is the coupled fermion-scalar system, then 

pE^em^n^pn, pn°^XF and XF°^a~^, (8) 

and we obtain 

^ln(e) ^Ixf de 
3 dln{xF) 3 e dxF 
This is the basic equation for computing w in our model, and it gives 



leFKoxl-3i2-y){l-y) 
3 eFKoxl + {2-y){l-y) 



: \ . I ' (10) 
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FIGURE 1. w vs. log(jCf ) for 8 values of Kq. 



It follows immediately that w > — 1 in our model. 

From equations (4,10), we can compute w as a function of xp, given ^o- Figure 1 
shows some results using a log scale for Xf. 

Note that w approaches close to — 1 as the density decreases (remember the direction 
of decreasing density is the direction of expansion of the universe) and then departs 
sharply towards zero. At large Xf, it is clear that w approaches +1/3 as it should for a 
relativistic gas of fermions. The value goes to zero at zero density. 

For large values of the minimum value of w is very close to — 1 and is reached near 
the value Xf = xf,i = {3/Ko) , which can be obtained from an analytic approximation. 
More precisely, and empirically jc^m ^ (3.82/^o)^^^- 



TABLE 1. The minimum value of w, the 
point Xf,,n at which this is achieved, and the 
approximate location Xf i 



log A'o 


log 


lOgX/7,1 


w 


6 


-1.809 


-1.841 


-0.94709 


9 


-2.806 


-2.841 


-0.99453 


12 


-3.806 


-3.841 


-0.99945 



The minimum value of w approaches —1 more closely as increases, and the low 
density recovery of w to zero becomes steeper as increases. 

It is clear from these results that we have a model in which the value of w is a function 
of the scale parameter a or the red-shift z, i.e. that we have an epoch dependent dark 
energy. But these results are in terms of dimensionless parameters and variables. To 
connect to the real world we must link these dimensionless numbers to the dimensionful 
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numbers which characterize it. As was shown in ref. ([3]), there are only a few, weak 
constraints on the actual parameter values. Moreover, very large values of Kq are possible 
even for very small values of if the range of the scalar is very large, corresponding to 
very small values of m^. Even if long-ranged, such weak interactions between fermions, 
especially neutrinos or those outside of the Standard Model altogether (such as the LSP) 
are exceptionally difficult to constrain by any laboratory experiments. 



CONNECTING WITH DARK ENERGY, AND FIXING 

PARAMETERS 

The energy density of dark energy is quoted as Pde = (3.20±0.4) x lO^'^^Gev'^ [6] If we 
define poE = tn\, then niE = 2.4 meV. It is important to remember that the conventional 
value of the dark energy density is derived on the assumption that it is constant during 
the evolution of the universe, and this is not the case in our model. Pending a study of 
the development of the universe using our model we use the present estimate of and 
apply it to the region w — 1 to estimate the relevant parameters in our model. 



,(0) 



PE^mrePn= ' " (U) 

To use this equation to estimate rn"n^ , first note that near w ^ —1, for large Kq, xp is 
small and we can use the appropriate approximations to get 

where, in the last equation, we use xp ~ xp.i at the minimum. Kq determines the 
relationship between m^^ and mp, ni^n'^ = [Ik^Kq^ nip. We choose Kq , taking care 
that the implied values of and nis are reasonable. For = 10^^ ^^n^ = 150 meV, 

in the neutrino range, and for Kq — 10^^, m"n^ — 150 GeV, in the range expected for 
neutralinos. What are the implications for the other parameters? One can immediately 
estimate the density of the neutral fermions as 47 x 10^ cm""^ and 3.5 x 10~^ cm~^ for 
these two cases. 

Because we are assuming that the neutral fermions form a homogeneous back- 
ground, it is appropriate to set the range of the scalar field to be the scale param- 
eter of the universe at the relevant time, which we will take to be z = 1, i.e. nis ~ 
(7 X 10^ lightyears) ~ 3 x lO^-'*^ meV. With this value of nis we can obtain the im- 
plied value of g^/{An): 3 x 10^^^ and 3 x 10^^^ in these two cases. Even the largest 
coupling is far too weak to be constrained by terrestrial experiments. 

An alternative way to proceed would be to assume a present density of the neutral 
fermions, po, and thus a present value of the Fermi momentum, kpQ, at a scale parameter 
flo> use akp = ao^FO and aq = (1 + z)a. Equation (1 1) then gives 

[i+zy Po 
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To give an explicit example, assume that the present density of neutral fermions is 
100 cm^, characteristic of neutrinos, and extrapolate this to a density of 800 cm~^ 

at z = 1. This gives a value ml^^ = 9.2eV, and Xf = 3 x 10^^, = 9.7 x 10^, and the 
even smaller coupling constant, g^/ {4-n) = 9 x 10^^^. 

At this value of the minimum value of w is —0.8, at the l.5o level from the value 
— 1.05 ±0.18 of the ESSENCE supernova survey[13], so even these parameters are not 
excluded. 

It is impossible to make further progress without finding additional ways to constrain 
the parameters. Given their extreme values the most promising approach is to use 
additional theoretical input, which we leave as a challenge for future work, by us and 
others. 



RESULTS — z DEPENDENCE OF w 

Now that we have some parameters we can convert the determination of w as a function 
of the dimensionless parameter Xf to a dependence of w on the red shift z. To illustrate 
the results we have selected the values 10^ and 10^ for and set the minimum of w to 
occur at z = 1 . 

These results are illustrated in figure 2. 




FIGURE 2. w vs. z for A^o ^ 10^ Kq = 10** and Kq ^ lO'^. 

Note that there is a slow variation of w with z at epochs earlier than z^^i^i at which the 
minimum occurs, but a rapid change for z between the present epoch (z = 0) and that of 
the minimum. 

With the present uncertainty in data we can choose to place the value of Zj^j^ at any 
value of z between, say z = 0.5 and z = 1.2 without being in conflict with the present 



Epoch Dependent Dark Energy 



August 10, 2009 



6 



data, and have chosen to show the results at Zm\n = 1 as illustrative results. 



CONCLUSIONS 

In an ideal future world one could imagine that w{z) is indeed observed to have a 
minimum value at some Zmin, and the the value of w-^jj^, close to but greater than 
— 1 is accurately known. The fact that the minimum value of w is near —1 shows that 
^0 > 10^'^, and the actual value of w^^, if known accurately enough, would determine 
the value of ^o- Knowing the value of allows us to predict a value of ps at the 
minimum w, which provides a good test on the model. Then a selection of parameters 
can place that minimum near an appropriate z value, z^[^ . 

Our prediction of a characteristic z dependence of w and pE leaves us vulnerable to 
developments in the precision of measurements and the extension of observations to 
larger and (with more difficultly) to smaller z, and we await with interest the results of 
the proposed experiments, such as those described in ref. [6]. 
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